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SUMMARY 


The objective of this work is to develop an analytical method that can provide insight 
into the delamination mechanism and to investigate the fatigue delamination 
characteristics. The approach is to derive a finite-element method, apply the method to 
layups designed to test the method and to innovative layups designed for improved 
properties, and then validate by specimen tests. 

A simple quasi-three-dimensional finite-element model has been developed to study the 
free-edge delamination in composite laminates subjected to uniaxial tension, bending, 
and torsion. This model was used to study laminates with different Poisson's ratios 
between plies and different stacking sequence. It was found that for a laminate under 
bending, minimizing the Poisson's ratio difference between plies or stacking the primary 
load-carrying plies far from the mid-plane can significantly reduce the interlaminar 
stresses and strain energy release rates. Unlike the tension case, the total strain energy 
-release rate (Gt) and the Mode I strain energy release rate (Gj) for laminates subjected 
to bending decrease as delamination length increases, reaching a lower bound. 

A family of quasi-isotropic graphite/epoxy laminates was investigated to determine the 
strain energy release rate at various interfaces in the laminate. In addition, a 90° 
graphite ply was replaced with a glass ply to create a hybrid laminate in the quasi- 
isotropic layup. A comparison of the strain energy release rates of all-graphite and 
hybrid layups indicated that the hybrid laminates show a significant reduction in the 
total strain energy release rate at all the interfaces and in the Mode I strain energy 
release rate at the 0/90 interface of the various layups. 

For the hybrid coupons, a 20% increase in the static onset-of-delamination strength and 
a 10% increase in ultimate strength resulted from replacing a 90° ply of graphite with a 
glass ply in the [±30 2 /90 3 /+30 2 ] T graphite laminate. These results confirmed the author's 
previous work. The test results of all-graphite laminates were included to serve as a 
baseline for comparison in both static and fatigue tests. For a fatigue life greater than 
106 cycles, the maximum fatigue load for this type of coupon cannot exceed 59% of the 
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static onset-of-delamination strength for the baseline coupon and 51% of same 
delamination strength for the hybrid laminate. It was also found that the hybrid 
laminate lowers the fatigue delamination onset. 

Coupons of S2/SP250 glass/epoxy laminates were tested to establish the delamination 
onset for Mode I and mixed-mode edge-delamination fatigue fracture. From the Mode I 
static edge-delamination test, the onset-of-delamination strength was used in the 
equations to determine interlaminar fracture toughness with and without consideration 
of the thermal effects of curing. It was found that including the residual stress due to 
curing increased the interlaminar fracture toughness 125%. The equation developed 
by Whitney et al. was used to calculate the fracture toughness that included the curing 
effect. Since the onset of delamination cannot be directly observed in this type of 
coupon, a delamination length 3% as long as the laminate width was selected as the 
criterion for the onset of delamination during fatigue. This standard was then used for 
all the coupons tested in this study and is recommended as superior to a loss-of-stiffness 
.standard. 

In both the Mode I and mixed-mode glass/epoxy laminate, the relationship between load 
and onset of fatigue delamination was established. The test results indicate that to 

achieve 1 x 106 cycles of fatigue life, the maximum fatigue load cannot be applied 
above 60% of the static onset-of-delamination load for the Mode I laminates and 43% 

for the mixed-mode laminates. 


xiv 



FOREWORD 


This report describes the development of a finite element model that is capable of 
identifying the delamination characteristics of laminates under various loadings. 
Fatigue delamination behavior is also depicted. This work was conducted by Bell 
Helicopter Textron, Inc. (Bell) between January 1986 and June 1988 in fulfillment of 
National Aeronautics and Space Administration Contract NAS1-18199. The contract 
effort was financially supported by the Army Aerostructures Directorate at the Langley 
Research Center. Dr. T. K. O'Brien was technical monitor. 

The project engineer of this program was Dr. Wen S. Chan of the Research Structures 
Group at Bell. He was assisted in the development of the finite element model by Dr. 
Ozden O. Ochoa during her leave of absence from Texas A&M University. The coupon 
test program was performed in the Methods and Materials Laboratory at Bell. 
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INTRODUCTION 


Composite laminate materials have played an important role in the overall improvement 
of airframe and rotor system structures. One of the most important, yet least 
understood, factors affecting the strength and durability of composites is interlaminar 
stress. The flapping capability of the composite rotor, for example, is limited by 
delamination resulting from the interlaminar stresses at the free edges and ply drops. 
Improvement in the free-edge delamination strength of composites depends on a 
fundamental understanding of the delamination characteristics of the structural 
laminates. Therefore, development of an analytical method that can provide insight 
into the delamination mechanism under complex loadings is of foremost importance to 
the analyst. 

Early efforts to understand the free-edge delamination problems focused on 
development of an analytical method to determine interlaminar stresses and strain 
energy release rates and on investigation into the delamination characteristics of 
laminates under extension. Pipes and Pagano (1970) solved the elasticity problem by a 
finite-difference method. Wang and Crossman (1977), using three-noded triangular 
elements, developed a quasi three-dimensional finite-element approach for modeling the 
cross sections of laminates. Raju and Crews (1981) took the same approach but used 
eight-noded isoparametric elements instead. Pagano and Soni (1983) used a global-local 
model to obtain the interlaminar stress components by lumping plies far from the 
interface of interest with an equipollent system. Wang and Choi (1982) used 
Lekhnitskii's stress function and an eigenfunction expansion to obtain a solution for the 
interlaminar stresses in a laminate. Kassapoglou and Lagace (1985) used the force 
balance method and the principle of minimum complementary energy to obtain the 
interlaminar stresses by assuming the shapes of the stresses. Using lamination theory, 
Sandhu and Sendeckyj (1986) calculated the delaminating moment coefficient as the - 
measure of delaminating tendency of laminates. All the works cited above addressed 
laminates subjected to extension only. 

Salamon (1978) extended Pipes and Pagano's finite-difference method to laminates under 
pure bending. Armanios and Rehfield (1986) included transverse shear deformation 
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through rotation, grouped plies together away from the interface of interest, and 
applied the ply elasticity equation to obtain the stresses between plies. Both tension 
and bending cases were studied. Their model provided a reasonable estimate for the 
interlaminar shear components but not for the normal stress component, since no 
thickness strain was included in the model. No work, to the author's knowledge, has 
been done on a laminate under torsion. 

The strain energy release rate has been widely adopted by analysts because of its 
success in characterizing delamination. Pagano and Pipes (1973) observed the 
interlaminar normal stress governing the delamination strength. O'Brien (1980) devised 
the laminate coupon test to investigate the onset and growth of delamination in 
graphite/epoxy laminates under uniaxial extension. He developed a simple expression 
for computing the total strain energy release rate, Gf The components of the total 
strain energy release rate, Gi, Gu, and Gm, associated with opening, shear, and tearing 
failure modes, respectively, were also based on a quasi three-dimensional finite-element 
analysis and the crack-closure technique (Rybicki and Kanninen, 1977). Armanios and 
Rehfield (1986) divided the delaminated region of the laminate and estimated G?, Gu, 
and Gm by the crack-closure method. Then Gi was obtained by subtracting Gu and Gm 
from Gj. Only uniaxial extension and bending loads were considered in their studies. 

This report covers both the analytical and experimental work performed in a study of 
the delamination characteristics of laminated composites. In the analytical work, an 
integrated two-dimensional finite-element model capable of incorporating various loads, 
including axial, bending, and torsion, is developed to calculate the interlaminar stresses 
as well as strain energy release rates in laminates. The study addresses the formulation, 
implementation, and verification of the model. In order to understand delamination 
characteristics, parametric studies of the effect on interlaminar stress and strain 
energy release rate of Poisson's mismatch between layers, stacking ' sequence, 
delamination length, and delamination opening for various laminates due to bending are 
included. Preliminary results on laminates subjected to torsion are also included. In the 
experimental work, fatigue delamination characteristics of graphite/epoxy, glass/epoxy, 
and their hybrid laminates are investigated. The work is focused on the establishment 
of fatigue thresholds for Mode I and mixed-mode edge delamination for glass and 
graphite/glass hybrid laminates. 
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INTEGRATED FINITE-ELEMENT ANALYSIS FOR FREE-EDGE PROBLEM 


A conventional two-dimensional finite-element model is not capable of analyzing the 
free-edge stress problem, and a full three-dimensional finite-element model is not 
efficient. Therefore, a quasi three-dimensional finite-element model capable of 
incorporating various loads, including tension, bending, and torsion, has been developed 
and is described in this section. 

FINITE-ELEMENT FORMULATION 

Displacement Function 

An isoparametric quadrilateral plate element with eight nodes and three degrees of 
freedom per node is used to model the cross section of the laminate in the y-z plane (see 
Figure 1). The displacement field is defined with the customary parabolic interpolation 
functions Nj, as given below: 


u(y,z) = £ N.U. 


v(y,z)= £ N.V. 


w(y,z) = £ N.W. 


Ni in the element coordinates (r,s) are expressed as 

Ni = i(l + r)(l + s)-i(l-r 2 )(l + s)-i(l-s 2 )(l + r) 

N 2 = i(l-r 2 )(l + s) 

N 3 = i (1-r) (1 + s) - i (1-r 2 ) (1 + s) - * (1— s 2 ) (1-r) 

N 4 = i(l-s 2 )(l-r) (2) 

N 5 = i(l-r)(l-s)-i(l-s 2 )(l-r)-i(l-r 2 )(l-s) 

N 6 = i(l-r 2 )(l-s) 

N 7 = i (1 + r) (1-s) - i (1-r 2 ) (1-s) - i (1-s 2 ) (1 + r) 

Ng = i(l-s 2 )(l + r) 
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ONLY SHADED AREAS ARE 
CONSIDERED IN THE FINITE- 
ELEMENT MODEL. 



Figure 1. Schematic of a quasi three-dimensional finite-element model. 




The strain field based on the displacement functions of Eq. (1) takes the following form: 




^z 

*xz 
* 


r 


yz 


*y 

V J 




3u/dx 

dv/dy 

dw/dz 

dv/dz + dw/dy 
d\l/dz + dw/dx 
du/dy + dv/dx 

V y 


[B] 

6X24 


{u} = 

24X1 


{^}o 


(3) 


where the [B] matrix represents the partial derivatives of the interpolation functions 
with respect to coordinates and {u} is the total number of nodal displacements, that is, 


Vi 

Wi 

• 

{u} = / . \ (4) 

• 

U 8 

V 8 

l w «J 

In order to successfully incorporate "bending" loads into the planar element, the 
displacement field of Eq. (1) is altered to the following expressions: 

U(x,y,z) = € 0 x + kxz + u(y,z) 

V(x,y,z) = Cxz + -v(y,z) (5) 

W(x,y,z) = - ^ kx 2 - Cxy + w(y,z) 


where € 0 is the uniform extension, and k and C are the longitudinal and twisting 
curvatures, respectively, as obtained from 
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( 6 ) 


K 


X 


aV 


ax‘ 


= k 


a 2 w 



a 2 w 

K = - 2 = 2C 

*y axay 


Thus, a constant bending curvature through the laminate thickness is enforced without 
bending-extension coupling. It is also noteworthy to look into the case of applying pure 
twisting curvature, which represents the pure torsion load. For this case K xy (i.e., C) 
will have a prescribed value where € 0 and k of Eq. (5) can be expressed as 


{€} = {€}o + {€}l 


(7) 


with {€} 0 being the strain vector of Eq. (3) and {€}l being defined as 


{€}l 


^kz + ^ 
0 
0 
0 

-Cy 

Cz 






( 8 ) 


It is important to reference two other displacement fields proposed for the bending 
problem, namely, those of Armanios and Rehfield (1986), and Ye and Yang (1986). The 
approach of Ye and Yang, formulated independently of the present authors' work, 
modifies the displacement field of the same eight-noded isoparametric element to 
incorporate longitudinal curvature only. The induced twist of laminates that bending 
contributed to the displacement was totally neglected. Therefore, their model can only 
be applied to the [0/90] type of laminates. The displacement field suggested by Rehfield 
is used in a closed-form solution rather than in a finite-element solution and includes a 
shear deformation through the rotation 0(y) (see Table 1). Comparisons of the three 
different displacement formulations and the corresponding strains are presented in 
Tables 1 and 2, respectively. It is immediately apparent that the present formulation 
extends itself to torsion problems, in addition to bending problems. The presence of 
zero strains, i.e., € z = € xy = € xz = 0, in the Rehfield approach is not desirable for 
accurate presentation of deformation. 
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TABLE 1. COMPARISON OP DISPLACEMENT FIELD FORMULATIONS 


Displacements 

Chan & Ochoa 

Ye & Yang 

Rehfield 

U 

u(y,z) + kxz 
+ € 0 x 

u(y,z) + kxy 

kxz 

V 

v(y,z) + Cxz 

v(y,z) 

v(y) + z0(y) 

w 

w(y,z)-}kx 2 

-Cxy 

w(y,z) - £kx 2 

w(y) - ^kx 2 


TABLE 2. 

COMPARISON OF STRAINS CORRESPONDING TO DISPLACEMENT 
FIELD FORMULATIONS 

Strains 

Chan & Ochoa 

Ye & Yang 

Rehfield 

€* 

€ 0 + kz 

kz 

kz 


dv 

dv 

8y 

3v zd(3 

ay” + dy 

^Z 

3w 

~dz 

0W 

dz 

0 

<*y 

— + Cz 
3y + Uz 

du 

dy 

0 

^XZ 

dll n 

I? - ^ 

d\l 

dz 

0 

^yz 

dv 3w 

~dz + dy 

du aw 

dz + 8y 

. dw 

P + Hy 
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Equilibrium Statement 


The present finite-element model formulation is based on the minimization of total 
potential energy. The strain energy of the system (all assembled elements) is expressed 
as 

0 = - £ [ {o} T {€}dVol (9) 

2 e = l Jvo1 

where N is the total number of elements and the volume integral is over an element or 
with the use of the constitutive relations 


N 

0 = z 


e = l 


~ {€ } T [ D] {€} dVol 

vol 2 


Substituting Eq. (7) into Eq. (10) yields 


« = I f \ (<«>J + <«£) ®] ({<>„ + («} L )dVoi 

e = 1 


( 10 ) 


(ii) 


Upon rearranging terms, 


Q 


N 

- I 

e- 1 


vol 


\ (<«ftDl <«)•„ + {€>r [D] {€} l 


+ {€}[ [D] {€} & + {€}* [D] (€} l ) dVol (12) 


Now, recall that from Eq. (3), {6} 0 can be represented by 


{e} o = [B]{u} 


(13) 
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Substituting Eq. (13) into Eq. (12) and collecting terms yields 


Si = 


I [ l ({u} t [B] t [D][B]{u} + {u} T [B] T [DJ{C} 

e= 1 J Vol 2 V 


+ {€}J [D][B]{u} + {e£[D]{<E} L ) dVol ( 14 ) 


On the other hand, the potential energy due to the general applied loads is 


N 

W 2 | - {u} T {f} dVol 

e= 1 


(15) 


Thus, the total potential energy is obtained by the summation of Eqs. (14) and (15). The 
minimization of the total potential energy with respect to the nodal displacements 
(generalized coordinates) provides the equilibrium statement for the system, that is, 

an _ a(Q+n>) _ ( 16 ) 

du S u 

The resulting equilibrium statement, in terms of the stiffness matrix and the load 
vector, is presented below: 


y [ [B] T [D][B] dVol){u} + (y f [B] T [D]{£} dVol 

e»i Vol ) y- Jvol L 


= i | i 

e = 1 J Vol 


(f}dVol 




V 




[K] 

Stiffness matrix 
for the assembled 
system 


(17) 



{f} 

Contribution from Acting on the 

the loading terms system in general 


The second term on the left side of Eq. (17) is obtained from the contributions of the 
curvatures included in the displacement field of Eq. (5). 


9 



Then, the final form of the equilibrium statement is given as 


I [ [B] T [D][B]dVol){u} = £ f ({f}-[B] T [D]{e} L ) dVol (18) 

e= l 1 Vol / e = i J VolV L / 


or, in short, 


[K] {u} = {F*} (19) 

As can be observed, the effect of curvature can be interpreted as a "bending load" 
effect. Therefore, one does not need to re-evaluate the stiffness matrix for the 
isoparametric element with an extensional displacement field, as given in Eq. (1). The 
bending effect can simply be incorporated by redefining the load vector, as shown in 
Eqs. (18) and (19). 

Representation of Loading Terms 

From classical lamination theory, the moments can be expressed as 


M x = Bii€ 0 + Djjk + D 12 Ky + 2DieC 

(20) 

My = Bi2^o + Di2k + D22Ky + 2D26C 

(21) 

M xy = BieC 0 + D^k + D 26 Ky + 2Dg6C 

(22) 


where € 0 is elongation in the x-direction, and k, K y , and C are the curvatures about x, y, 
and the twisting axes, respectively. Positive sign convention for k is concave downward. 

For a problem of pure bending about the x axis, M y = M xy = 0. Therefore, simultaneous 
solutions of Eqs. (21) and (22) yield the expression for 2C (twisting curvature) as a 
function of k (longitudinal curvature) and € 0 (elongation): 
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2C - _ n n 2 [*o( B 12 D 26“ B 16 D 22) + k ( D 12 D 26 D 16°22) 


(23) 


22 66 26 


For symmetric laminates under pure bending, the above expression reduces to 


2C 


D 22 D 66 °26 


k ( D 12 D 26 - D l6 D K ) 


(24) 


Then, the strain field obtained from the bending contributions, {€}l, can be expressed as 
a function of k, only, by substituting Eq. (24) into (8). 


A similar discussion is valid for pure torsion load about the x axis with M x = M y = 0. 
Then, for a prescribed value of C, kand K y are given as 


k = — 2C 


D 16 D 22- D 12 D 26 \ 
D 22 D 12- D ?2 ' 


(25) 



D 11 D 26 D 12 D 16 
D 22 D 11 - °?2 


(26) 


It should be noted that the bending stiffness, Dy, used here was obtained from 
lamination theory for laminates with no delamination. To account for laminates with 
edge delamination, the bending stiffness should be re-evaluated. 


Boundary Conditions 

Because of symmetries inherent in the laminates under consideration and in the loading 
condition across the cross section, the region to be modeled is dependent on the loading 
case. The three cases, shown in Figure 1, are modeled as follows: 
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Uniaxial loading A portion of the quarter cross section close to the free 

edge, where interlaminar stress exists, is modeled with 
symmetric boundary conditions imposed along the y and z 
axes. 

Bending The total through-the-thickness direction is modeled with 

symmetric boundary conditions imposed about the z axis. 

Torsion The entire cross section is modeled, creating antisymmetric 

displacements and stresses. 

The origin of the modeled region is constrained against any displacement, i.e., u(0,0) = 
v(0,0) = w(0,0) = 0. 

VERIFICATION OF BENDING CAPABILITY 

A thick laminate, [0/902/0]i6s> 1 in^ wide and 0.64 in. thick, was modeled with 120 
elements. The material was T300/5208, with the properties shown in Table 3. The z 
axis was taken as an axis of symmetry; thus, only a half-width of the laminate was 
modeled. The applied loading was in the form of prescribed curvature K x = k = 0.1. The 
outermost layers were 0.02 in. thick, which is 6.25% of the total laminate thickness. 
This geometry and the loading conditions ensured the existence of uniform stress-strain 
fields in the outermost thin layers. To assess the results obtained from the computer 
program for the bending load, another model was used to obtain the stress field for a 
thin [0/90]s laminate, 0.02 in. thick (same thickness as the outermost layer of the thick 
laminate), subjected to uniaxial tension. Figures 2 and 3 show, respectively, the normal 
and shear stresses obtained from the bending model and the uniaxial tension model at 
the upper 0/90 interface. The good agreement between the two models lends confidence 
in the programming of the bending subroutine developed for the computer code. 

FINITE-ELEMENT RESULTS 

In this section, the finite-element code based on the formulation described in the 
previous section is verified. The computer code is then used to investigate the effect of 
interlaminar stresses resulting from a Poisson's ratio mismatch between plies and the 
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stacking sequences in laminates, and the effect on strain energy release rate of changes 
in the magnitude of delamination openings and delamination sizes. A comparison of the 
present model with the approximate closed-form model is presented. The effect of 
strain energy release rate in hybrid laminates is also discussed. 

To determine the strain energy release rate, the bending stiffness used was calculated 
from the laminate with no delamination, since the length of the delamination considered 
in this study was small compared to the laminate width. 

The material constants for the laminates used in this study are listed in Table 3. 


TABLES. MATERIAL CONSTANTS 


Parameter 

IM6G/3501-6 

S2/SP250 

S2/CE9000 

T300/5208 

Ei 

24.8 Msi 

7.3 Msi 

8.09 Msi 

20.0 Msi 

E 2 

1.41 Msi 

2.1 Msi 

2.28 Msi 

2.1 Msi 

G 12 

0.9 Msi 

0.88 Msi 

1.06 Msi 

0.85 Msi 

V 12 

0.329 

0.275 

0.288 

0.21 

tply 

0.0074 in. 

0.0085 in. 

0.0074 in. 

0.005 in. 


Parametric Study 

The investigation of interlaminar stresses resulting from different Poisson's ratios 
between plies and from different stacking sequences is described below. 

Effect of Poisson's Ratio Mismatch 

The difference in Poisson's ratio mismatch between plies has a significant effect on 
interlaminar stresses of laminates subjected to tension (Chan et al., 1986). To study this 
effect on interlaminar stresses of laminates subjected to bending, a series of T300/5208 
laminates with a [02/©23s layup were used. The 0 was varied from 30° to 90° in incre- 
ments of 15°. The applied loading was in the form of prescribed curvature K x = k = 0.1. 
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(ksi) 


{0/90 2 /0] 16S 



Figure 2. Variation of interlaminar normal stress o zz with y at 0/90 interface 
under bending load. 



figure variation oi interlaminar snear stress o yz with y at u/90 interlace 
under bending load. 
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Figure 4 shows the distribution of normal stress, a zz , across the width at the upper 0/0 
interface of the [O 2 /02]s laminates where the 0° plies exhibit positive longitudinal stress 
(<j x ). As shown in the figure, the laminates with 0 = 90° and 75° give rise to an 
interlaminar normal tensile stress at the free edge. However, the laminates with 0 = 
15°, 30°, 45°, and 60° have an interlaminar normal compressive stress at the edges. 
The tensile and compressive interlaminar normal stresses at the free edge can be 
explained by the difference in Poisson's ratio between the plies. 

The variation in Poisson's ratio with respect to the fiber orientation of the material 
considered in this study is depicted in Figure 5. For the 0 = 90° or 75° laminate, there is 
less deformation in the transverse direction than in its neighboring 0° ply, so the edges 
of the laminate at this interface tend to curve away from the midplane of the laminate. 
As a result, interlaminar normal tensile stress is induced. For 0 = 60° or less, the edges 
of the laminate at the upper 0/0 interface curve inward, since the less deformed 0° ply 
opposes the effort of its neighboring ply to curve outward, and the result is compressive 
-stress in the thickness direction. It is noted that at the lower 0/0 interface, all the 
stresses will be reversed, since the fiber stress of the 0° plies is under compression. 

Figure 6 shows the distribution of interlaminar normal stress, a zz , through the thickness 
near the free edge. The figure confirms the zero a zz at the outer surfaces and at the 
midplane, where no longitudinal stress exists because of bending. It is also shown that 
the stress behaves, as expected, as an antisym metrical distribution. 

The interlaminar shear distribution, a yz , along the upper 0/0 interface for the 0 = 30°, 
60°, and 90° laminates is shown in Figure 7. In this figure, the maximum peak for 
interlaminar shear at the 0/0 interface occurs in the 0 = 30° laminate, which is the angle 
where the maximum difference in Poisson's ratio occurs. 

The interlaminar shear variation through the thickness near the free edge, illustrated in 
Figure 8, indicates that the peak values occur at the interface and at the midplane 
locations. 
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200 


O 15° 



- Figure 4. Variation of interlaminar normal stress a zz with y at z/H = 2 for [O 2 /02Js 
laminates under bending load. 



0 (deg) 


Figure 5- Variation of Poisson's ratio v*y with 0 in T300/5208 graphite/epoxy laminate 


o 30° 



z/H 


Figure 6. Variation of interlaminar normal stress a M through-the- thickness for [O2/02]s 
laminates at (b - y)/H = 1 under bending load. 



Figure 7. Variation of interlaminar shear stress o yz with y at z/H = 2 for [Oo/Bols 
laminates under sending load. 



Figure 8. Variation of interlaminar shear stress o yz through-the-thickness for [O2/02ls 
laminates at (b - y)/H = 1 under bending load. 

Effect of Stacking Sequence 

Two families of laminates, (09/±452)s S2/SP250 glass/epoxy and (±45/0/90)s IM6/3501-6 
graphite/epoxy, were used to study the influence of stacking sequence on interlaminar 
stresses and strain energy release rates. The material constants for these materials are 
given in Table 3. 

Laminates of S2/SP2S0 glass/epoxy with three different stacking sequences were 
fabricated: [09/±452ls> [07/±45/0/±45/0]s, and [05/±45/03/±45/0] s . The ply thickness, H, 
was 0.0085 inch. The laminates were subjected to a constant curvature of k = 0.1 and 
plots made of the resulting interlaminar stresses. These plots, shown in Figures 9 and 
10, depict the distribution of normal and shear components, respectively, through the 
thickness at one ply distance away from the free edge. It is clear that all the peaks for 
normal stress occur at the first ±45 ply interface from the outer surface of the 
laminate. However, the peaks for shear are located at the first 0/45 interface from the 
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Stacking Sequence 
|0^±4S 2 | s 



Figure 9. Effect of stacking sequence on interlaminar normal stress in laminates subjected to bending. 



Stacking Sequence 
[0$/ ± 45}]$ 



s 


Figure 10. Effect of stacking sequence on interlaminar shear stress in laminates subjected to bending, 


outer surface. The results also strongly indicate that, of the three laminates, [09/±452]s 
has the smallest peak values for both normal and shear. The calculated strain energy 
release rates for these three laminates with a delamination at the first 0/45 interface 
from the bottom surface are shown in Table 4. The results imply that Gj and Gn, as well 
as Gt> can be significantly reduced if all the 0° plies are placed near the outer surface. 

TABLE 4. EFFECT OF STACKING SEQUENCE ON STRAIN ENERGY 
RELEASE RATES OF S2/SP250 LAMINATES® 


Strain Energy Release Rate (in- lb/in 2 x 10' 3 ) 


Laminate 

Gi 

Gn 

Gill 

Gt 

[O 9 / ± 45 2 ] S 

3.033 

3.471 

0.817 

7.321 

[ 0 7 / ± 45/0/ ± 45/OJs 

14.170 

15.123 

1.518 

30.811 

[ 05 /± 45 /O 3 / ±45/0] S 

24.1.96 

26.528 

1.358 

52.082 


*k = 0.1, a/H = 8. Delamination at first 0/45 interface from bottom. 


The stacking sequences for the IM6/3501-6 graphite/epoxy laminates were [0/±45/90]s, 
[±45/0/90]s, and [90/±45/0]s. The stress fields are displayed in Figures 11 and 12. 

Figure 11 shows a high peak for a zz at the (45/-45) interface (z/H = -2) in the [90/±45/0]s 
laminate. However, a high peak for the shear, cr yz , shown in Figure 12, is located at the 
90/45 interface (z/H = ±3) of the [90/±45/0]s laminate. 

Comparison of Present Results 

The [0 8 /90 4 ] s laminate of T300/5208 graphite/epoxy studied by Armanios and Rehfield 
(April 1986) was selected for comparison with the present finite-element results. Since 
their theory can only provide interlaminar shear, interlaminar normal stress was not 
considered in this study. 

Figures 13 and 14 show the interlaminar shear (o yz ) distribution along the 0/90 interface 
from the top of the laminate when subjected to bending and tension, respectively. It is 
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Figure 11. Effect of stacking sequence on interlaminar normal stress o zz through-the- 
thickness at (b-y)/H = 1 for IM6/3501-6 laminates under bending load. 

(90/ ♦ 45/-45/0]s 

[ + 457*45/0/90]$ 



Figure 12. Effect of stacking sequence on interlaminar shear stress a yz through-the- 
thickness at (b-y)/H = 1 for IM6/3501-6 laminates under bending load. 


99 



Figure 13. Comparison of results for interlaminar shear stress o yz at z/H = 4 for the 
[0g/90 4 ]s laminates under bending load. 

T300/5208 Gr/gp Laminate with Creek 



Figure 14. Comparison of results for interlaminar shear stress o yz at z/H = 4 for the 
[0 8 /90 4 ] s laminates under tension load. 


23 


noted that the interlaminar shear stress zone depicted by both methods is comparable 
for laminates subjected to tension. However, the finite-element method predicts a 
larger zone than Rehfield’s method in the case of bending load. This is because the 
contribution of the shear strain, in the z direction to the interlaminar shear, a yz , is 
not considered in Rehfield's solution. 

Both methods were used to obtain the strain energy release rates of the same laminate. 
Figure 15 shows that the finite-element results are comparable to Rehfield’s results. In 
both cases the G calculations were based on a delamination having a length of six ply 
thicknesses located at the 0/90 interface. 



Gj G| G|| 

Figure 15. Comparison of present model and Rehfield's results for strain energy release 
rates in [0 8 /90.i]s T300/5208 laminates under bending load. 

Height of Vertical Delamination Opening vs Strain Energy Release Rates 

A [90 4 /0/±45]s IM6/3501-6 graphite/epoxy laminate was subjected to a constant 
curvature to study the relationship between strain energy release rate and the height of 
a vertical (z-direction) delamination opening at the 0/45 interface. Both the opening 
and the length of the delamination were predetermined in the analysis. Three different 
ooening heights were considered. The results, shown in Table 5, indicate that the height 


of the delamination opening has little if any effect on the G calculation. It was 
therefore concluded that the use of a mathematical crack with zero height to simulate 
the actual delamination in the analysis would not result in error. 


TABLE 5. EFFECT OF VERTICAL CRACK OPENING HEIGHT ON 
[90 4 /0/±451s IM6/3501-6 LAMINATE WITH K x = 0.1 


Crack 

Opening 

Grp 

(in-lb)/in 2 

Gi 

(in-lb)/in 2 

Gii 

(in-lb)/in 2 

Gm 

(in-lb)/in 2 

0.0 H* 

12.0x10-3 

4.12x10-3 

3.22 x 10-3 

4.87 x 10-3 

0.02 H 

12.2x10-3 

4.12x10-3 

3.19x10-3 

4.91 x 10-3 

0.1 H 

12.2x10-3 

4.12x10-3 

3.22x10-3 

4.89 x 10-3 


*Ply thickness 


Delamination Length vs Strain Energy Release Rates 

To determine how the strain energy release rates change as the length of a delamination 
increases in the y-direction, a [90 4 /0/±45]s laminate of IM6/3501-6 graphite/epoxy was 
subjected to a constant curvature of k = 0.1. Since the laminate was found to have the 
highest interlaminar stresses at the lower 0/45 interface prior to delamination, a 
delamination was assumed to occur at that interface. The results of the calculation are 
shown in Figure 16. For very short cracks, i.e., a/H = 1 to 5, the behavior of the 
individual G's does not follow a consistent pattern. However, at a delamination length 
greater than eight times the ply thickness, the strain energy release rates, Gi, Gn> Gm> 
and G-r* level out and become constant in value. This behavior is different from that 
observed in the tension case. For tension, the strain energy release rate increases as 
the delamination length increases, reaching a plateau at a certain length. To better 
understand the G's behavior, the forces and displacements used for G calculation are 
shown in Figure 17. It is noted that F z and U z decrease significantly as delamination 
length increases. 
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Figure 16. Change in strain energy release rates with increasing delamination length for 
a [90 4 /0/±45]s IM6/3501-6 graphite/epoxy laminate. 



16 

-14 
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Figure 17. Forces and displacements used in calculation of strain energy release rates 
for [90 4 /0/±45]s IM6/3501-6 graphite/epoxy laminate. 
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Effect of Hybrid Laminate 


A family of quasi-isotropic laminates with a baseline of [45/-45/0/90]s and off-axis 
intervals of 15° was examined. The total strain energy release rate and the Gi, Gn, and 
Gni components were evaluated at three interfaces. For three of these laminates, the 
90° graphite ply was replaced with a glass ply of S2/CE9000, creating a hybrid laminate; 
these laminates were [45/-45/0/90]g, [0/90/-45/45]s, and [-45/45/90/0]s. As can be 
observed in Figures 18 through 20, the hybrid laminate consistently demonstrated a 
reduction in the total strain energy release rate. It is also noted that the [-45/45/90/0]s 
layup (Figure 20) has the highest strain energy release rate of the three layups. 
Moreover, at the 0/90 interface of this layup, the hybrid laminates exhibit a significant 
reduction in the Mode I strain energy release rate, Gi. Since Gi is considered to be a 
major contributor to static delamination, hybridization of soft 90°-ply laminates can 
significantly increase the delamination strength. Figures 21 through 24 depict the strain 
energy release rate components for the rest of the off-axis quasi-isotropic family. Of 
these, the [-45/45/90/0] layup has the highest Gi at the 45/90 interface. 



N H N H N H 

Figure 18. Comparison of strain energy release rates at various interfaces for the 
[±45/0/90]s laminates, with and without the hybridizing 90° ply. 
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N- baseline 90° ply 
H - hybrid 90° ply 


V/S/S/SS/A 




Figure 19. Comparison of strain energy release rates at various interfaces for the 
[0/90/+45]s laminates, with and without the hybridizing 90° ply. 


N- baseline 90° ply 
H - hybrid 90° ply 


XrssssssssA 


Figure 20. Comparison of strain energy release rates at various interfaces for the 


[+45/90/0] s laminates, with and without the hybridizing 90 ply. 








- 15/75 - 60/- 15 30/-60 


Figure 21. Variation of strain energy release rates at various interfaces for the 
[30/-60/-15/75]s laminates. 
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- 30/60 - 75/-30 15/-75 


Figure 22. Variation of strain energy release rates at various interfaces for the 
[15/-75/-30/60]s laminates. 
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Figure 23. Variation of strain energy release rates at various interfaces for the 
[-15/75/-60/3Q]s laminates. 
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Figure 24. Variation of strain energy release rates at various interfaces for the 
[-30/60/-75/15]s laminates. 


Preliminary Results on Torsion 


Three layups of 16-ply IM6/3501-6 graphite/epoxy laminates - [± 45 ) 45 , [±452/0 2 /902]s> and 
[0]gs - were studied under pure torsion loading. Two delaminations were assumed to be 
present, each extending in from the edge for a length equivalent to eight ply 
thicknesses. One was located at the top interface on one side of the laminate and the 
other was at the bottom interface on the opposite side. Their strain energy release 
rates for laminates subjected to a 0.1 twisting curvature are given in Table 6. As 
shown, a pure Mode III was achieved in the 0°-ply lamina, and the ±45 laminate yielded 
the highest Gt. The data also indicate that as the percentage of ±45 plies increases, Gn, 
Gin, and Gt increase. Further studies on delamination characteristics have also been 
performed by Chan and Ochoa (1987). 


TABLE 6. STRAIN ENERGY RELEASE RATES FOR LAMINATES 
SUBJECTED TO A 0.1 TWISTING CURVATURE 


G (in-lb/in 2 ) 

[±45] 4S 

[±452/0 2 /90 2 ]s 

[0]gs 

G, 

0.1307 

0.3332 

0 

Gn 

1.0066 

0.0181 

0 

Gin 

1.9259 

0.6872 

0.2806 

Gt 

3.0632 

1.0389 

0.2806 
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EXPERIMENTS 


The development of materials with delamination resistance has been of interest to both 
material manufacturers and designers for some time. While manufacturers have been 
devoted to improving toughness by altering the properties of the material, designers 
have focused on improving toughness by using conventional materials in innovative 
laminate constructions. New materials, such as toughened thermoset and thermoplastic 
matrices, have shown a significant increase in resistance to static delamination, but 
little improvement in resistance to fatigue (O'Brien, March 1987). On the other hand, 
innovative laminate constructions, using conventional materials and devices such as an 
inner layer (Chan et al., 1986), critical ply termination near the free edge (Chan and 
Ochoa, 1988), and a hybridized soft 90° ply, have exhibited a even greater improvement 
in static delamination resistance * far beyond the capability of the current toughened 
matrices. Both inner layer and critical ply termination concepts have also successfully 
demonstrated a remarkable increase in fatigue resistance. The hybrid laminates have 
not yet been tested for fatigue delamination resistance. 

The purpose of this experimental program was to investigate the fatigue delamination 
characteristics of graphite/epoxy, glass/epoxy, and their hybrid laminates. The specific 
objectives were to establish the fatigue delamination thresholds for the Mode I and 
mixed-mode fatigue fractures of glass and graphite/glass hybrid laminates. 

The static and fatigue thresholds for delamination were determined for both Mode I and 
mixed-mode fractures to characterize the material behavior of S2/SP250 glass/epoxy 
laminates. The fatigue characteristics of a hybrid laminate consisting of an AS4/3501-6 
graphite/epoxy laminate with a 90° ply of S2/CE9000-9 glass/epoxy were investigated. 
Both 3501-6 and CE9000-9 matrix are 350°F cure resins. 


STATIC AND FATIGUE TEST MATRIX 

Edge delamination tension tests were first introduced by O'Brien (1980) to investigate 
the characteristics of delamination. Using layups such as [±302/90 3 /?302 Jt> O'Brien was 
able to measure delamination before final failure by reducing the amount of strain 
necessary to initiate the delamination. In this layup, the delamination grew at both the 
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top and bottom -30/90 interfaces, producing a mixed-mode fracture. Whitney and 
Knight (1985) modified this layup by embedding strips of a homogeneous material in the 
midplane of the laminate to force the delamination to grow at that plane, thus 
producing the Mode I fracture. 

The test matrix and laminate layups, as well as the materials used in each test type, are 
outlined in Table 7. It should be noted that the hybrid coupons contain AS4/3501-6 
graphite/epoxy plies for ±30° laminae and S2/CE9000-9 glass/epoxy plies for 90° 
laminae. A Tinius Olsen test machine was used for the static tests and an MTS 880 
servohydraulic machine for the fatigue tests. The static tests were conducted to obtain 
the load at which static onset of delamination takes place, since this load is used in the 
determination of the maximum fatigue loads. 


TABLE 7. COUPON TEST MATRIX 


Test Type 

Coupon Layup 

Material 

Load Type 

Test 

No. of 
Tests 

Mode I 

[30/-30 2 /30/90] s 

S2/SP250 

Static 


2 




Fatigue 

f = 1Hz 

9 




R = 0.1 

f = 5 Hz 

9 

Mixed-Mode 

( ± 3 O 2 / 9 O 3 / T 302)t 

S2/SP250 

Static 


2 




Fatigue 

f = 1Hz 

9 




R = 0.1 

f = 5 Hz 

9 

Mixed-Mode 

[ ± 30 2 /902Gl/* 3021 t 

AS4/3501-6 

Static 


2 


Hybrid 

S2/CE9000-9 

Fatigue 

f = 1Hz 

9 




R = 0.1 

f = 5 Hz 

9 


COUPON PREPARATION AND GEOMETRY 

Panels measuring 12 inches on each side were fabricated from each laminate. To 
simulate initial delaminations, 0.55-inch-wide strips of half-mil thick Kapton were 
placed in some panels. These strips were located 1 inch apart in the 0° direction 
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at the midplane of the laminate. Each panel was vacuum-bagged and autoclave-cured, 
with a final cure temperature of 250°F for glass laminates and 350°F for the graphite 
hybrid laminate, since SP250 is a 250°F curing matrix and both CE9000-9 and 3501-6 are 
350°F curing matrices. The panels with the Kapton strips were cut along the centerline 
of each strip so that each coupon had a Kapton strip approximately 0.25 inch wide along 
its free edges. Before and after cutting, the panels and coupons were checked by C- 
scan for any imperfection that might have occurred during fabrication and cutting. 

The coupon configurations for the laminates are shown in Figure 25. All coupons were 
approximately 1.5 inches wide and 10 inches long. Glass-laminate end tabs were 
attached to all coupons. 

TEST PROCEDURE 

Coupons with the Kapt on-simulated initial delamination at both of the free edges were 
•tested for Mode I fracture, while coupons without the Kapton were tested for mixed- 
mode fracture. 

Static Tests 

In the static tests, the coupon response during loading was monitored by means of a 1- 
inch MTS extensometer mounted on two V-notched aluminum plates attached to the test 
coupons 1 inch apart. An X-Y plotter was used to record the stress-strain behavior 
during loading. In the Mode I static test, two C-ring gages developed by Boeing Aircraft 
Corporation of Seattle were used to monitor the onset of delamination. In addition, a 
flashlight was shone through the back side of the coupons to enable visual detection of 
any delamination that might have occurred. In the mixed-mode test, checks for 
delamination were made by periodic visual inspections through a microscope. A typical 
coupon with the gages installed is shown in Figure 26. 
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Figure 26. Setup for static tests. 


Fatigue Tests 

The fatigue tests were conducted at 1 and 5 Hz under maximum fatigue loads of 90%, 

85%, 75%, and 60% of the static onset-of-delamination strength and a stress ratio of R 
(=Omin/omax) = 0#1 * The cou P ons were inspected for delamination growth at every 500 
cycles from 0 to 5,000 cycles, at every 1000 cycles from 5,000 to 10,000 cycles, and at 
every 5,000 cycles from 10,000 to 100,000 cycles. If a delamination occurred, the 
inspection was conducted only at 5,000, 10,000, and 1,000,000 fatigue cycles. Dye 
penetrant was applied at the edge of the coupon to enhance detection of the 
delamination area. For glass/epoxy laminates, ink was used; for graphite/epoxy 
laminates, zinc iodide, which is opaque to x-rays, was used. Photographs were taken at 
the time of the inspection to document the delamination size. For graphite/epoxy 
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laminates, x-radiographs were taken. During all the tests, coupon edges were 
periodically checked by touch to see if any rise in temperature could be detected. Such 
a rise in temperature could indicate the release of energy and, hence, a change 
occurring within the coupon. Although subjective in nature, these tactile inspections 
are useful in initiating closer inspections of an objective nature. 

TEST RESULTS AND DISCUSSION 

The test results for each type will be discussed in the following section. The fracture 
toughness values (strain energy release rates) will also be presented. 

Mode I Edge-Delamination Tension Test 

Static Test Results 

•The static test results are listed in Table 8. The onset of delamination was measured 
from the stress-strain curve. Although the C-ring gage described in the previous section 
was mounted on the coupon to monitor the onset of delamination, there was some 
difficulty in measuring the onset value because of slippage of the gage. Therefore, the 
onset of delamination was measured at the change in slope on the stress-strain curve. A 
typical load-strain curve is plotted in Figure 27. The first slope change is caused by 
transverse cracks, and the subsequent slope change is caused by delamination. The 
measured value at the onset of delamination was verified by visual inspection (light 
shining through the back side of the coupon). 

The unidirectional material constants listed in Table 9 were used to calculate the 
fracture toughness. The stress-free temperature used for estimating the thermal effect 
during curing was assumed to be 182°F. Subtracting the room temperature (72°F) 
results in a AT of -110°, which was the value used in calculating the Mode I fracture 
toughness (Gic)* The measured static onset-of-delamination strain and the calculated 
fracture toughness (see Appendix), with and without the thermal effects of curing, are 
included in Table 8. The results indicate that Gic is significantly affected by the curing 
process. The initial delamination (Kapton) located at the midplane of the coupon divides 
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TABLE 8. STATIC G IC FOR S2/SP250 GLASS/EPOXY 


Calculated 

Measured Interlaminar Fracture Toughness 

Onset of Delamination (in-lb/in2) 


Coupon 

UltimateStrength 

(ksi) 

Strength 

(ksi) 

Strain 

(pe) 

With 

Thermal Effects 

Without 
Thermal Effects 

i 

45.23 

18.79 

6300 

0.630 

0.286 

2 

45.07 

17.33 

6000 

0.587 

0.260 

3 

41.00 

16.01 

5700 

0.545 

0.234 

Avg 

43.77 

17.38 

6000 

0.587 

0.260 



Figure 27. Typical load-strain curve for Mode I edge-delamination coupon. 
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TABLE 9. UNIAXIAL MATERIAL PROPERTIES 
FOR S2/SP250 GLASS/EPOXY LAMINATE 


Material 

Ei 

(Msi) 

mm 


Via 

Q 1 

(pe/°F) 

a 2 

(pe/°F) 

tply 

(in) 

S2/SP250 

6.31 

2.50 

0.60 

0.25 

10.73 

42.12 

0.0095 

S2/CE9000-9 

8.09 

2.28 

1.06 

0.288 

- 

- 

0.0074 

AS4/3501-6 

19.3 

1.62 

1.02 

0.288 

5 

25 

0.0052 


the symmetric laminate into two unsymmetric sublaminates in the delamination region. 
Hence, the curing effect will play an important role in this type of test. 

Fatigue Test Results 

Figures 28, 29, and 30 show delamination growth, Aa, beyond the initial delamination 
(Kapton strip) for coupons subjected to 90%, 85%, and 75% of the static onset-of- 
.delamination load, respectively. The. delamination growth, Aa, was measured from the 
maximum extent at the same six selected locations for all the coupons. The selected 
locations were located approximately at the middle length of the coupon, and 2 inches 
up and down from the middle length on each side of the coupon. This procedure was 
then used for all the coupons tested in this study. The growth is plotted as a function of 
number of fatigue cycles and has been normalized by coupon width, b. The delamination 
growth was measured during the inspection intervals. The figures also show data from 
the two vibration rates used, 1 Hz and 5 Hz. Although the indication is not strong, 
there appears to be a tendency, under all three loads, for the delamination to grow 
faster at 1 Hz than at 5 Hz. 

The fatigue delamination onset was selected as that point at which the delamination 
growth beyond the initial length reached 3% of the coupon width. This value was 
determined by a least-squares curve fit of the data shown in Figures 28, 29, and 30. The 
delamination growth shown in these figures is the growth beyond the width of the 
Kapton strip. Table 10 shows the resulting fatigue delamination onset for the Mode I 
coupons. In Figure 31, these values are plotted against the applied load, which is given 
as a fraction (or percentage) of the static onset-of-delamination, a/o 0 . A straight line, 
least-squares fitted on a semi-logarithmic scale, is drawn through the data and 
extrapolated to 1 x 106 cycles, which was assumed to be the endurance life. To achieve 
this, the load, a/a 0 , for Mode I would be approximately 0.60. 
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TABLE 10. FATIGUE DELAMINATION ONSET FOR MODE I EDGE-DELAMINATION 

GLASS/EPOXY COUPONS 


Coupon 

Width X Thickness 
(in X in) 

* o 
°0 

Freq 

(Hz) 

Fatigue Delamination 
Onset 
(cycles) 

Avg 

(cycles) 

1 

1.507 X 0.098 

0.90 

1 

8,793 


2 

1.502 X 0.096 

0.90 

1 

7,338 

6,871 

3 

1.507 X 0.100 

0.90 

1 

4,483 


4 

1.494 X 0.099 

0.90 

5 

8,702 


5 

1.493 X 0.097 

0.90 

5 

10,240 

8,791 

6 

1.496 X 0.096 

0.90 

5 

7,432 


7 

1.494 X 0.098 

0.85 

1 

15,860 


8 

1.495 X 0.100 

0.85 

1 

9,671 

15,067 

9 

1.504 X 0.100 

0.85 

1 

19,670 


10 

1.503 X 0.098 

0.85 

5 

34,580 


11 

1.491 X 0.099 

0.85 

5 

9,047 

19,082 

12 

1.500 X 0.098 

0.85 

5 

13,620 


- 13 

1.496 X 0.101 

0.75 

1 

72,390 


14 

1.501 X 0.098 

0.75 

1 

45,530 

>72,640 

15 

1.505 X 0.099 

0.75 

1 

>1 X 105 


16 

1.500 X 0.100 

0.75 

5 

>1 X 105 


17 

1.500 X 0.100 

0.75 

5 

63,550 

>87,850 

18 

1.502 X 0.101 

0.75 

5 

> 1 X 105 



a o/o„ = Maximum fatigue load/static onset of delamination; 
static delamination strength = 17.38 ksi 
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Figure 31. Effect of load level on fatigue life of Mode I edge-delamination 
glass/epoxy coupons. 
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Mixed-Mode Edge-Delamination Test for Glass/Epoxy Coupons 


Static Test Results 


Table 11 lists the static test results for the glass/epoxy mixed-mode coupons. The 
onset-of-delamination strength shown in the table will serve as a guide in determining 
the fatigue load. The values for interlaminar fracture toughness, G c , were calculated by 
the following equation: 




E 


lam 


- E 


) 


where 


e c = axial strain at the onset of delamination 
t = laminate thickness 

Eiam = laminate axial stiffness calculated from laminated plate theory 
E* = axial stiffness of a fully delaminated laminate 


For the [±302/903 /t302]t laminate, delamination occurred at the interface between the 
-30° and 90° plies. E* was calculated as follows: 


E = 


8E (±30)s + 3E 90 
11 


A detailed discussion of the above method can be found in O'Brien et al. (1982) and 
Standard Tests for Toughened Resin Composites (1983). 


TABLE 11. STATIC TEST RESULTS FOR MIXED-MODE 
EDGE-DELAMINATION COUPONS 


Coupon 

Onset of Delamination 
Strength Strain 

(ksi) (pe) 

Interlaminar Fracture 
Toughness, G c 
(in-lb/in 2 ) 

i 

21.28 

7600 

0.818 

2 

21.55 

7300 

0.754 

Avg 

21.42 

7450 

0.786 
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Fatigue Test Results 


As described in the Mode I test, the fatigue delamination onset was determined at the 
life when the delamination reached 3% of the lamination width. Table 12 presents the 
fatigue delamination onset for the S2/SP250 glass/epoxy coupons subjected to mixed- 
mode fracture. In Figure 32, the fatigue delamination onset values are plotted against 
the applied load, o/a 0 , and a least-squares-fitted straight line drawn through the data 
and extrapolated to 1 x 106 cycles. This shows that for a delamination life greater than 
10® cycles, the maximum fatigue load cannot be above 43% of the static onset-of- 
delamination strength. Figure 33 shows, for comparison, one of the coupons after 5,000 
cycles and again after 100,000 cycles. 

Mixed-Mode Edge-Delamination Test for Hybrid Coupons 

In this section, the results of static and fatigue tests on both graphite and graphite/glass 
.hybrid coupons are reported. The all-graphite/epoxy results generated under a Bell 
IR&D program serve as a baseline for comparison with the hybrid coupon results. 

Static Test Results 

Table 13 summarizes the static test results of coupons with and without a glass hybrid- 
ization. The data indicate a 20% increase in delamination strength and a 10% increase 
in ultimate strength for the hybrid laminate compared to the baseline. This improve- 
ment was observed by Chan and Wang in 1986. The improvement in delamination 
strength can be explained by the reduction in strain energy release rate shown by the 
analysis. For both the baseline and hybrid laminates, delamination was observed as wan- 
dering at the interfaces between the -30° and 90° layers and between the 90° layers. 

Fatigue Test Results 

Tables 14 and 15 summarize the fatigue test results for the graphite and graphite/glass 
hybrid coupons, respectively. It is evident that the hybrid coupons have 
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TABLE 12. FATIGUE DELAMINATION ONSET FOR MIXED-MODE EDGE- 
DELAMINATION GLASS/EPOXY COUPONS 


Coupon 

Width X Thickness 
(in X in) 

*0 

°0 

Freq 

(Hz) 

Fatigue Delamination 
Onset 
(cycles) 

Avg 

2 

1.486 X 0.105 

0.90 

1 

500 


3 

1.498 X 0.106 

0.90 

1 

572 

536 

4 

1.502 X 0.105 

0.90 

5 

463 


5 

1.495 X 0.106 

0.90 

5 

303 


6 

1.505 X 0.104 

0.90 

5 

702 

489 

8 

1.498 X 0.103 

0.75 

1 

8,201 


9 

1.481 X 0.105 

0.75 

1 

5,469 

6,840 

10 

1.500 X 0.104 

0.75 

5 

4,579 


11 

1.485 X 0.106 

0.75 

5 

3,438 

4,009 

12 

1.492 X 0.104 

0.70 

1 

9,647 


13 

1.499 X 0.104 

0.70 

1 

8,993 

9,320 

14 

1.498 X 0.104 

0.70 

5 

10,680 


15 

1.486 X 0.105 

0.70 

5 

8,579 

9,630 

16 

1.496 X 0.106 

0.60 

1 

62,280 


17 

1.496 X 0.105 

0.60 

1 

80,220 

71,250 

18 

1.486 X 0.105 

0.60 

5 

56,780 

56,780 


a o/o 0 = Maximum fatigue load/static onset of delamination; 
static delamination strength = 21.42 ksi 


a lower threshold!) life than the baseline coupons. Figures 34 and 35 show the fatigue 
delamination onset plotted against appplied load, o/o 0 , and a least-squares-fitted straight 
line drawn through the data and extrapolated to 1 x 106 cycles. For an endurance life 
greater than 10« cycles, the maximum fatigue load cannot exceed 59% of static onset- 
of-delamination strength for the baseline coupon and 51% of same delamination strength 
for the hybrid laminate. It is concluded that although the hybrid laminate can improve 
both the static delamination and ultimate strength, it lowers the fatigue delamination 
onset 
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TABLE 13. MIXED-MODE DELAMINATION TENSION TEST FOE GRAPHITE/ 
EPOXY WITH AND WITHOUT A GLASS HYBRIDIZATION 


Coupons 

Dimension 
Width X Thickness 
(in X in) 

Onset of Delamination 
Strength Strain 
(ksi) (pe) 

Ultimate Strength 
(ksi) 

All-graphite 

1.498 X 0.061 

26.05 

3800 

64.25 


1.493 X 0.059 

30.08 

3900 

63.57 

Avg 


28.07 

3850 

63.91 

Hybrid 

1.503 X 0.060 

33.05 

5200 

68.93 


1.502 X 0.061 

34.11 

4840 

71.16 

Avg 


33.58 

5020 

70.05 


TABLE 14. 

FATIGUE DELAMINATION ONSET FOR BASELINE EDGE-DELAMINATION 

GRAPHITE COUPONS 

Coupon 

Width X Thickness 
(in X in) 

a o 

°o 

Freq 

(Hz) 

Fatigue Delamination 
Onset 
(cycles) 

Avg 

i 

1.494 X 0.058 

0.90 

1 

276 

276 

2 

1.492 X 0.058 

0.90 

5 

78 

78 

3 

1.493 x 0.059 

0.75 

1 

7,458 


4 

1.495 X 0.058 

0.75 

1 

4,480 

5,969 

5 

1.492 X 0.060 

0.75 

5 

27,140 


6 

1.494 X 0.058 

0.75 

5 

17,810 

22,475 

7 

1.493 X 0.061 

0.70 

1 

23,050 


8 

1.493 X 0.060 

0.70 

1 

>1 X 105 

>61,525 

9 

1.495 X 0.059 

0.70 

5 

23,640 


10 

1.493 X 0.063 

0.70 

5 

>1 X 105 

>61,820 

11 

1.492 X 0.062 

0.60 

1 

> 1 X 105 

>100,000 

12 

1.493 X 0.063 

0.60 

5 

> 1 X 105 

> 100,000 


a o/o 0 = Maximum fatigue load/static onset of delamination; 
static delamination strength = 28.07 ksi 
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TABLE 15. FATIGUE DELAMINATION ONSET FOR HYBRID EDGE-DELAMINATION 

GRAPHITE/GLASS COUPONS 


Coupon 

Width X Thickness 
(in X in) 

*0 

°0 

Freq 

(Hz) 

Fatigue Delamination 
Threshold 
(cycles) 

Avg 

1 

1.499 

X 

0.061 

0.90 

1 

102 

102 

2 

1.500 

X 

0.060 

0.90 

5 

93 

93 

3 

1.495 

X 

0.060 

0.75 

1 

4,058 


4 

1.498 

X 

0.064 

0.75 

1 

2,620 

3,190 

5 

1.500 

X 

0.060 

0.75 

1 

2,892 


6 

1.499 

X 

0.060 

0.75 

5 

3,570 


7 

1.497 

X 

0.060 

0.75 

5 

2,261 

3,291 

8 

1.500 

X 

0.059 

0.75 

5 

4,041 


9 

1.500 

X 

0.062 

0.70 

1 

3,693 


10 

1.499 

X 

0.060 

0.70 

1 

11,770 

10,208 

11 

1.499 

X 

0.059 

0.70 

1 

15,160 


12 

1.498 

X 

0.060 

0.70 

5 

14,770 


13 

1.501 

X 

0.060 

0.70 

5 

18,790 

14,657 

14 

1.501 

X 

0.060 

0.70 

5 

10,410 


15 

1.501 

X 

0.059 

0.60 

1 

24,640 


16 

1.501 

X 

0.060 

0.60 

1 

38,580 

31,610 

17 

1.501 

X 

0.060 

0.60 

5 

243,400 


18 

1.494 

X 

0.064 

0.60 

5 

>100,000 

>100,000 


a o/o 0 = Maximum fatigue load/static onset of delamination; 
static delamination strength for baseline laminate = 28.07 


Figure 36 presents x-radiographs of two pairs of baseline and hybrid coupons after 
100,000 cycles under a load of 90% on the static onset of delamination. One pair was 
subjected to a rate of 1 Hz, and the other pair to a rate of 5 Hz. As can be seen, there 
are no significant differences in the delamination growth at this load level. 

In contrast to the coupons under 90% load, Figure 37 shows that delamination grows 
faster in the hybrid than in the baseline when a load level of 75% is applied at a rate of 
1 Hz. The x-radiographs in this figure were also taken after 100,000 cycles. 
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Figure 36. Comparison of 1- and 5-Hz baseline and hybrid graphite coupons 
after 100,000 cycles under 90% load. 










CONCLUSIONS 


A simple quasi-three-dimensional finite-element model has been developed and 
successfully used to study the free edge stress as well as the edge delamination 
characteristics of laminates subjected to uniaxial, bending, and torsional loads. This 
analysis and specimen tests of conventional as well as innovative layups have led to the 
following conclusions: 

1. For a laminate under bending, minimizing the Poisson's ratio difference 
between plies or stacking the primary load-carrying plies far from the mid- 
plane can significantly reduce the interlaminar stresses and strain energy 
release rates. 

2. Unlike the tension case, the total (Gt) and Mode I strain energy release rates 
for laminates under bending decrease as delamination length increases and 
reaches a lower bound. 

3. Gn and Gm, as well as Gt, increase as the percentage of 45° plies increases 
for a laminate subjected to torsion. 

4. The static onset-of-delamination strength and strain, as well as the ultimate 
strength, are higher in the graphite/epoxy laminates with the hybridized 90° 
glass ply than in the baseline all-graphite laminates. This improvement can 
be attributed to the reduction in strain energy release rate in the hybrid 
laminate, which is shown in the analysis. However, this improvement under 
static loading is not duplicated in fatigue test. The fatigue delamination 
onset is lower in the hybrid laminate than in the all-graphite laminate. 

5. For S2/SP250 glass/epoxy laminates, the calculated Mode I interlaminar 
fracture toughness (Gic) for the Mode I edge-delamination test is up to 125% 
higher when the thermal effects of curing are taken into account. 
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RECOMMENDATIONS 


The following recommendations are made for future research: 

1. Extension of the present work to increase understanding of the failure 
mechanism and delamination characteristics of edge delamination in 
unsymmetric laminates and multiple delaminations. 

2. Quantitative assessment of the effect of differences in Poisson's ratio and 
shear strain mismatch between two adjacent plies on delamination 
characteristics of laminates under tension, bending, and torsion loads. 

3. Investigation of the frequency effect on fatigue delamination threshold by 
testing at frequencies ranging from 0.1 to 30 Hz, to cover creep at both low 
and high frequencies. 

4. Development of an analytical model that includes the frequency effect to 
predict the fatigue delamination threshold. 


PRECEDING PAGE BLANK NOT FILMED 
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APPENDIX 


CALCULATION OF CRITICAL STRAIN ENERGY RELEASE RATE FOR 
THE MODE I EDGE-DELAMINATION TEST COUPON 


The calculation of strain energy release rate for the Mode I edge-delamination tension 
test was originally derived by Whitney (1985). In reviewing his work, errors were found 
in equation (14) of the original paper. In this section, details of the derivation are not 
included. Instead, only the equations used for the calculation are presented. 


The critical strain energy release rate with the residual thermal effect included can be 
written as follows: 

°,c=[N“.- T ) +T ,]K- T -n-[ E *«- T ) +T 2 ]( E « + 'W)| 


where 
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a = delamination length 
b = half of coupon width 
h = half of coupon thickness 
ti = ply thickness 

kt*> ply reduced stiffness 

W ij 


aj = ply thermal coefficients 

AT = temperature difference from stresss free condition 

If the thermal effect due to curing is neglected, Eq. (A-l) can be simplified to O'Brien's 
equation: 


o,c= K 2 ( e s -e'> 


62 



REFERENCES 


Armanios, E. A., and Rehfield, L. W., "Interlaminar Analysis of Laminated Composites 
Using a Subliminate Approach," presented in the 27th AIAA SDM Conference, May 1986. 

Armanios, E. A., and Rehfield, L. W., "Interlaminar Fracture Analysis of Composite 
Laminates Under Bending and Combined Bending and Extension," presented at the ASTM 
Eighth Symposium on Composite Materials Testing and Design, April 1986. 

Chan, W. S., Rogers, C., Cronkhite, J. D., and Martin, J., "Delamination Control of 
Composite Rotor Hubs," J. American Helicopter Society , August 1986. 

Chan, W. S., and Ochoa, 0. O., "Assessment of Free-Edge Delamination Due to Torsion," 
Proceedings of the American Society for Composites, 2nd Technical Conference, 
September 23-25, 1987. 

Chan, W. S., and Wang, A. S. D., "A Study on the Effect of the 90° Ply in Matrix Cracks 
in Composite Laminates," presented in the 27th AIAA SDM Conference, May 1986. 

Kassapoglou, C., and Lagace, P. A., An Efficient Method for the Calculation of 
Interlaminar Stresses in Composite Materials, Part 2 - Solution Using Assumed Stresses 
Shapes, TELAC Report 85-5, Technology Laboratory for Advanced Composites, Dept, of 
Aeronautics and Astronautics, Massachusetts Institute of Technology, March 1985. 

O’Brien, T. K., "Characterization of Delamination Onset and Growth in a Composite 
Laminate," Damage in Composite Materials, K. L. Reifsnider, ed., ASTM STP 775, 
American Society for Testing and Materials, 1982. 

O'Brien, T. K., Johnston, N. J., Morris, D. H., and Simonds, K. A., "A Simple Test for the 
Interlaminar Fracture Toughness of Composites," SAMPE Journal, Vol. 18, No. 4, 
July/August 1982. 


63 



O'Brien, T. K., Mixed- Mode Strain- Energy- Release Rate Effects on Edge Delamination 
of Composites , NASA Technical Memorandum 84592, NASA Langley Research Center, 
January 1983. 

O'Brien, T. K., "Delamination Durability of Composite Materials for Rotorcraft," 
presented at the 1987 NASA/Army Rotorcraft Technology Conference, NASA Ames 
Research Center, Moffett Field, California, March 17-19, 1987. 

Pagano, N. J., and Pipes, R. B., "Some Observations on the Interlaminar Strength of 
Composite Laminates," Int. J. of Mechanical Science, Vol. 15, 1973, pp. 679-688. 

Pagano, N. J., and Soni, S. R., "Global-Local Laminate Variational Model," Int. J. Solids 
Structures, Vol. 19, No. 3, 1983, pp. 207-228. 

Pipes, R. B., and Pagano, N. J., "Interlaminar Stresses in Composite Laminates Under 
-Uniform Axial Extension," J. Composite Materials , October 1970, pp. 538-548. 

Raju, I. S., and Crews, J. W., Jr., "Interlaminar Stress Singularities at a Straight Free 
Edge in Composite Laminates," Computers and Structures, Vol. 14, No. 1-2, 1981, pp. 
21-28. 


Rybicki, E. F., and Kanninen, M. F., "A Finite Element Calculation of Stress Intensity 
Factors by a Modified Crack Closure Integral," Engineering Fracture Mechanics, Vol. 9, 
1977, pp. 931-938. 

Salamon, N. J., "Interlaminar Stresses in a Layered Composite Laminate in Bending," 
Fibre Science and Technology, Vol. II, 1978, pp. 305-317. 

Standard Test for Toughened Resin Composites, NASA Reference Publication 1092, 
revised edition, compiled by ACEE Composites Project Office, NASA Langley Research 
Center, Hampton, Virginia, 1983. 


64 



Wang, A. S. D., and Crossman, F. W., "Some New Results of Edge Effect in Symmetric 
Composite Laminates," J. Composite Materials , Vol. 2, 1977, pp. 92-106. 

Wang, S. S., and Choi, I., "Boundary-Layer Effects in Composite Laminates, Part 2 - 
Free Edge Stress Solutions and Basic Characteristics," J. Applied Mechanics , Vol. 49, 
1982, pp. 549-560. 

Whitney, J. M., and Knight, M., "A Modified Free-Edge Delamination Specimen," 
Delamination and Debonding of Materials , ASTM STP 876, W. S. Johnson, ed., American 
Society for Testing and Materials, Philadelphia, 1985, pp. 298-314. 

Ye, L., and Yang, B. X., "Interlaminar Stress Analysis in Symmetric Composite 
Laminates Under Bending," Proceedings of the International Conference on 
Computational Mechanics , G. Yagawa and S. N. Atluri, eds., May 25-29, 1986. 


65 


Standard Bibliographic Page 


1. Report No. 

NASA CR- 181824 

2. Government Accession No. 

3. Recipient's Catalog No. j 

i 

4. Title and Subtitle 

An Analytical and Experimental Investigation of Edge Delamination 
in Laminates Subjected to Tension, Bending, and Torsion 

5. Report Date j 

March 1989 ; 

6. Performing Organization Code ! 

7. Authorts) 

Wen S. Chan 

8. Performing Organization Report No. 

699-099-265 

9. Performing Organization Name and Address 
Bell Helicopter Textron, Inc. 

P.O.Box 482 

Fort Worth, Texas 76101 

10. Work Unit No. 

505-63-01-05 

11. Contractor Grant No. 

NASl-18199 

12. Sponsoring Agency Name and Address 
Langley Research Center 
National Aeronautics and Space Administration 
Hampton, Virginia 23665 

13. Type of Report and Period Covered 

Final, Jan 1986 - June 1988 

14. Sponsoring Agency Code 

15. Supplementary Notes 

Langley Technical Monitor: T. K. O'Brien 

16. Abstract 

An integrated two-dimensional finite element has been developed to calculate interlaminar stresses and strain 
energy release rates for the study of delamination in composite laminates subjected to uniaxial tension, 
bending, and torsion loads. This work addresses the formulation, implementation, and verification of the 
model. Parametric studies were conducted on- the effect of Poisson’s ratio mismatch between plies and the 
stacking sequence on interlaminar stress, and on the effect of delamination opening height and delamination 
length, due to bending, on strain energy release rate for various laminates. A comparison of strain energy 
release rates in all-graphite and graphite/glass hybrid laminates is included. The preliminary results of 
laminates subjected to torsion are also included. Fatigue tension tests were conducted on Mode I and mixed- 
mode edge-delamination coupons to establish the relationship between fatigue load vs. onset of delamination 
cycle. The effect on the fatigue delamination onset of different frequencies (1 and 5 Hz) was investigated for 
glass, graphite, and their hybrid laminates. Although a 20% increase in the static onset-of-delamination 
strength and a 10% increase in ultimate strength resulted from hybridizing the all-graphite laminate with a 
90° glass ply, the fatigue onset is lower in the hybrid laminate than in the all-graphite laminate. 
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